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Abstract 

In several applications, such as WENO interpolation and reconstruction [Shu C.W.: SIAM Rev. 51 (2009) 82-126], 
we are interested in the analytical expression of the weight-functions which allow the representation of the approx- 
imating function on a given stencil (Chebyshev-system) as the weighted combination of the corresponding approxi- 
mating functions on substencils (Chebyshev-subsystems). We show that the weight-functions in such representations 
[Muhlbach G.: Num. Math. 31 (1978) 97-110] can be generated by a general recurrence relation based on the ex- 
istence of a 1 -level subdivision rule. As an example of application we apply this recurrence to the computation of 
the weight-functions for Lagrange interpolation [Carlini E., Ferretti R., Russo G: SIAM J. Sci. Comp. 27 (2005) 
1071-1091] for a general subdivision of the stencil {xi-m_, ■•• , x;+m + } of M+ 1 := M_ +M + + 1 distinct ordered points 
into K s + 1 < M := M_ + M+ > 1 (Neville) substencils {x,_ M _+^, • ■ • . Xi+M+-K,+k,} (k$ e {0, • ■ ■ , K s }) all containing the 
same number of M - K s + 1 points but each shifted by 1 cell with respect to its neighbour, and give a general proof 
for the conditions of positivity of the weight-functions (implying convexity of the combination), extending previous 
results obtained for particular stencils and subdvisions [Liu Y.Y., Shu C.W., Zhang M.R: Acta Math. Appl. Sinica 
25 (2009) 503-538]. Finally, we apply the recurrence relation to the representation by combination of substencils of 
derivatives of arbitrary order of the Lagrange interpolating polynomial. 
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1. Introduction 

The Neville- Aitken algorithm [ 1 , pp. 204-209] constructs the interpolating polynomial on { jc,-m_ , • • ■ , *;+m + }, by 
recursive combination of the interpolating polynomials on substencils, with weight-functions which are also polyno- 
mials of x 111, pp. 204-209]. Carlini et al. |2|], working on the Lagrange interpolating polynomial in the context of 
centered (central) WENO schemes!^, recognized the connexion between the Neville algorithm JTl pp. 207-208] and 
the determination of the optimal [3] weight-functions. 

Definition 1.1 (Stencil). Let 

, m ( M + eZ:M:=M_+M + >0 
be a set of M + 1 distinct ordered real points. □ 
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Definition 1.2 (Neville substencils). Let X,_m ;+m + be a stencil (Definition II. lb and assume M > 2 in (fl}. The 
K s + 1 > 1 substencils 

/ N 3 < M - 1 := M_ + M + - 1 > 1 

X/-M-+* 8 ,/+M t -«:»+*» •- w-M-+t 8 » - " »^/+M t -isr 8 +* 8 J \ £■ e {0 ■ • • A"} 

each of which contains M - K s + 1 > 2 distinct ordered points and which satisfy 

Xj-M_+k s ,i+M + -K i +k s - X,_M_,i+M + (2b) 

* s =0 

Xi-M_+k s + l,i+M + -K s +k s +\ = {Xi-M_+kJ+M+-K s +K \ {Xi-M_+K }) U {*i+M + -A' s +«:s + l I VA: S € {0, • ■ • , K & - 1 } (2c) 

correspond to the A" s -level subdivision of x,-_m_,/+m + to substencils of equal length, each obtained from its left neigh- 
bour by deleting the leftmost point and adding 1 point to the right ( f2cb . □ 

The optimal weight-functions 07 >X; _ M _ , +M+ ,Ar s ,/t s (.x) in WENO interpolation satisfy ami 

Pi,Xi-u-MM + f> = J] °~i^i-M-MM + ,k„A (x) P/,x i _„_ +w+M+ _^ +Ss (Jc; /) Vx e M (3a) 

a; 



a; 



fc=0 



where Pi,x^ M _ i+M+ (x;f) is the Lagrange interpolating polynomial [1, pp. 183-189] of the real function / : K — > K on 
the stencil X,_m_,!+m + (Definition ! l.U . The optimal weight-functions o"/ x f M i+M ,^,jt s (x) (O correspond to the weight- 
functions in Miihlbach's theorem [5, Theorem 2.1, p. 100], where they were expressed in terms of quotients of 
determinants of interpolation-error functions, directly obtained by the Cramer solution [6, Proposition 5.1.1, p. 72] 
of error-eliminating linear systems 17, (13), p. 8489]. Since WENO approaches are based on nonlinear, with respect 
to the function f(x), modifications of the optimal weight-functions (O, we are particularly interested in analytical 
explicit expressions of the weight-functions. 

Carlinietal. Q (3.6,4. 10), pp. 1074-1079] gave the expression of the polynomial weight-functions o"/ iXKr l) . +rjr _i iis (x) 
for the (A^ = r — l)-level subdivision (Definition ll.2l ) of x,-_( r _iy +r (Defini tion ll.lt . This result was also confirmed 
by Liu et al. (2.2), p. 506] who further gave the analytical expression |0, (2.18), p. 511] for the polynomial 
weight-functions 07x ; .-iX*) f° r tne C^s = r)-level subdivision (Definition 1 1.21 ) of X,_, v+f (Definition ll.lt . In both 
cases it was observed ||2|,|4J] that Vx 6 [x,_i, x, + i] the linear weight-functions are positive (> 0), so that, because of the 
consistency relation Obi , the combination d3at of substencils is convex Vx e [x,_i, x, + i]. 

The purpose of the present note is to prove (Lemma [2. It a general recurrence relation [8, (70), p. 299] for weight- 
functions of an arbitrary family of functions, for which the (A" s = l)-level subdivision (Definition 1 1 ,2b is possible. 
As an example of application we apply this relation to the Lagrange interpolating polynomial (Proposition l3.lt . for 
an arbitrary level of subdivision (Definition 1 1.2t of a general stencil (Definition ll.lt . The explicit expression of the 
weight-funcions developed in Proposition 13. II is used (Proposition 13. 2t to study the convexity of representation 0J. 
Then, we apply the general recurrence relation (Lemma I2.lt to determine the weight-functions for the representation 
of the n-derivative of the Lagrange interpolating polynomial by combination of substencils (Proposition l4.2t . 



2. General recurrence relation for weight-functions 

We start by considering a more abstract case, where a general family of functions Pm-,m+( x ) depending on 2 
integer indices M ± € Z : M_ + M+ > 1 (which in practical applications may correspond to stencils; Definition ll.lt , 
are equipped with a 1 -level subdivion property, and develop a general recurrence relation for the weight-functions. 

Lemma 2.1 (Recursive generation of weight-functions). Let Pm-,m+( x ) be a family of real functions 

Pm.,m* '■ R -> K VM ± 6 z : M :- M- + M+ > 1 M+ > -M- (4a) 



2 



and assume that there exists an associated family of real weight-functions o~m_,m + xo( x ) (also defining o~m _,Af + ,l,l(x) := 
1 - o" m_,m + , i,o(*)) 



0"M_,Af+,l,oW + C"M_,M + , l,lW = 1 



VM ± e Z : M := M_ + M+ > 2 
Vfc s e{0,l} 

VM + eZ : M := M_ + M+ > 2 



defined everywhere in M except for a finite set of isolated points Sa- M M+ , C R, which may be empty, such that 



Vx e K \ S, 

VM ± e Z : M := M_ + M+ > 2 



Pm„,mA x ) — °~M-,M + , 1,0 W PAf_,M + -lW + CTM-.M+.l.lW Pm_-i,m+W 
r/ien the family of weight-functions defined recursively by 

min(X,-l,*s) 

I Vfc s e {(),■•■ ,2sT s } 



VM ± eZ : M_ + M + > 2 



(4b) 
(4c) 

(4d) 
(4e) 



4=max(0,jfc s -l) 



satisfies 



k,=0 



Vx e R \ (J (J 5, 

Z,=0 4=0 



PM-,M+( X ) — ^M-AAaW PM--k s ,M + -K ! ,+kX x ) 

VM ± e Z : M := M_ + M + > 2 
1 V^ s < M - 1 

Furthermore, for the values of{x, M ± , K s ] satisfying the conditions of (l4fb , 

^ Q"M-,m + ,a: s a(^) = 1 

jfcs=0 



(4f) 



(4g) 



Proof. Assume M ± e Z : M := M_ + M+ > 3 => (M_ - 4) + (M+ - 1 + € s ) = M - 1 > 2 W, e Z. Then, <gd]i 
applies to both functions /?m_-4,m+-i+4 (4 £ {0, 1}), and we have 



PM.-e s ,M + -i+eX x ) = Y\ 0"M_-4,M + -i+4,i,m s W P(M_-f s )-m s ,(*f + -i+4)-i+m s W j VM ± e Z i M_ + M + > 3 (5a) 

m s =0 ( 6 {0,1 1 

where cr M __f s s ,i , m s (x), being 1-level weight-functions, exist by ( I4"b"ll4dl >. Combining d4dll5al l. we have 



Vx e K \ 5o- M __ fs „ + _ 1+fs 



Ptt-JitJ< x ) ^ 0"M_,M + ,1/ S (x) C r Af_-4,M + -l+4,l,m s W pAf_-£ i -m S) Af + -2+4+m s (•*) 



4=0 
1 l 



\m*=0 



: ^ / , 0"Af_,M + ,l,4W 0"M_-4,M + -l+4,l,m s ( x ) PM_ 



-m s ,M + -2+4 +m 



(X) 



4=0 m s =0 
2 z' min(U s ) 



jt s =0U s =max(0,jt s -l) 

fM.,M t ,2,t,(-t) 



PM--k s ,M+-2+kX x ) 



(5b) 



Vx e R \ (5 £r „_, M+J U S aM _ M+ _ hi U ^„__ ljH+J ) VM ± e Z : M := M_ + M + > 3 => 2 < M - 1 
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which proves ( l4el l4fb . for K s = 2, because 

2-1 L, 1 

l^J \^_) 3°~M^-t s JU + -L s +t s ,\ ~ ( [^j ^0"M_-f s .M + +f s ,l) U ( l^j ^0-M_-f s ,M + -l+f s ,l) = (^M-.M+.l ^ ^0~M_,M+-U ^ ^f"M_ - 1 .*/+ , 1 ) (^ C ) 

L s =0 4=0 4=0 £ s =0 

To prove ( f4el [4fb VK S e {1, ■ • ■ , M - 1}, by induction, assume that ( f4el [4fb are valid for some K s - \ >2. Then 

K.-2 L. 



f s =0 



i s =o 4=0 (5d) 
{ VM ± e Z : M := M_ + M + >K s + \ 



with <x M M+ ^ s _] ti (x) in d5dl i defined by (l4el . Assuming ^T s < M — 1 in d5dV we can further subdvide each of the /f s 
functions /?m_-4,m + -c&: s -i)+4(- ;i in ( l5di i. once more, and we have by ( l4dl i 



f»_-CM t -K-l)+(,W = V 0"M_-4,M + -(if,-l)+A,l,m,W PM_-C s -m s ,M+-K^e,+m s ( x ) \ VM ± E Z : M_ + M + > # s + 1 

m,=o [ V€ s e {0,25T S - 1} 

(5e) 



where <TM_-€ s ,M + -(K s -i)+e s ,i,'n s (x), being 1 -level weight-functions, exist by ( l4bll4dl i. Combining ( T5dll5el i. we have 

K,-l ( 1 



,M+(x) ^ 0~M-,M+,K t -l,eS x ) ^ ^M_-4,M + -(tf s -l)+4,l,m s (*) f l Af_-4-m s ,Af + -X s +4+m 8 W 



4=0 

1 



\m s =0 



/ , / , 0"Af_,Af+,X s -l,4W 0"M--4,M + -(AT s -l)+4,l,m s (-^) P M--e,-m s ,M+-K s +i s +m,(x) 



4=0 m s =0 
A", / min(tf s -l,iy 



= ^ ^ o"m_,m + ,*: s -i,4W ^jf.-i.A-n-iWJA-f.W 

/t s =0V4=max(0,/t s -l) 

O" M_,M+,K s ,k s (x) 



PM_-k i ,M+-K a +k a ( x ) (5f) 



if. -2 L, 



Vx e R \ ( U U ^„__ fs ,„ + _ is+fsJ ) U ( |J S aM __ hM+ _ iKM ) VM ± e Z : M : = M_ + M + > JST S + 1 

' L„=0 4=0 4=0 ' 

This completes the proof of (14ft with weight-functions d4et . by induction. By d4eb . we have 



^ 0"M.,M*,K s ,kS. X ) ^ 0"M_,Af + ,X s -l,4W 0" J vf_-4,M + -(# s -l)+4,lA-4(-£) 



/f s min(A" s -l,fc) 



fe=0 



jt s =0 4=max(0,/t s -l) 
K s -\ 1 



(A.3)] V V ^ ^ 

- 2_i / i a "M-,M*£,-l.&W °"M_-4,M t -(jf s -l)+4,l,m s W 
4=0 m s =0 

= / \ o~m_,m+,k s -ia( x ) / j °"Af_-4,Af + -(^ s -l)+4,l,m a (-y) = / , O'tf-.M^-U.W (5g) 
4=0 V m s =0 J 4=0 

/e the sum of the weight-functions d4eb is the same VK S 6 {1, • • • , M — 1} (by induction). Since, by d4ct . ( 4g i holds for 
K s = 1, (|5gl suffices to prove d4gl V/f s €{!,••• ,M - 1}, by induction. □ 



3. Application to the Lagrange interpolating polynomial 

By Aitken's Lemma |H pp. 204-205], the Lagrange interpolating polynomial satisfies the 1-level subdivision 
property (l4bl |4cV with weight-functions which are linear polynomials, and therefore defined Vx € R, implying that 
Si,xi-M t+M ,k s — in ©. Application of Lemma |2~T1 to the Lagrange interpolating polynomial can be summarized in 
the following proposition. 



Proposition 3.1 (Weight-functions for the Lagrange interpolating polynomial). Assume the conditions of Definition ! 1.2 

Then, the weight-functions 07 x,_ M MM+ ,K„h( x ) ^ n ^ ne representation (01 of the Lagrange interpolating polynomial on 
X,-m_,/'+m + , Pi.xi-m- i+M+ (■*> f)> are rea l polynomials of degree K s with only real roots, expressed by 



,K„k, 



n 



(6a) 



M-J+M+ \X-i-M-+ks,i+M+-K s +ks 



V£ s e {(),•■■ ,tf s } VK S e {!,-•• ,AT- 1 := M_ +M + - 1) 



where the strictly positive real numbers 6/,x,_ m _,+m .a" s ,* s depend on the points of the stencil X,_m_./+m + (Definition ! 1.1b , 
anof are generated by the recurrence relation 



1 



X i+M + X i-M- 
min(A;-U s ) 



K s = 1 



4=max(0,i s -l) 

V£ s e {0,--- ,K S ] V£ s 6{1,--- ,M- 1 :=M_ +M + - 1} 
77ie weight-functions ( I6al l satisfy the consistency condition ( Bbl l ana? f/ie recurrence relation d4et . 
Proof. The case = 1 

1 



(6b) 



,lJt 



-(x - x, +m+ -<: s m) Vfc s e {0, 1 } 



(7a) 



X i+M + X i-M_ 

holds since it is exactly Aitken's Lemma [ 1, pp. 204-205]. Since © hold for K s = 1, VM ± eZ:M:=M_ + M + >2 
the family of Lagrange interpolating polynomials is equipped with the 1-level subdivison rule (l4al - l4cl i. and therefore 
satisfies the conditions of Lemma |2~T1 We can therefore apply d4eb to 07 x M _. +M ^^(jc). To obtain the simpler 
expressions (|6), assume that ( f6ab holds for some K s — 1 > 1, Then by Lemma lZTI 



0"/, Xi _ M _. i+M+ ,/(: s ,i: s (^) ^ ^ o"/,x,_ M _, i+M+ ,^ s -i,4(^) "/,x f _ M _ + f sJ+ „ + _ (Ks _ 1)+fs ,i,t s -f s ( J:: ) 



min(X s -l,ifc s ) 



£ s =max((U s -l) 
min(X s -l,ifc s ) 



2 ' ' 

f s =max(0,/fe s -l) ^ 

x (-1) 1 kl+ei Biju-H-n t j+M+-vc-w„iA-t, 



=(-1) 



]~ [ (X - X„) 

Xn^Xi-M--^e s j+M+-(K s -i)+es\^i-M-+e^+k s -e s j+M + -(K s -i)^e s -i+i^-e s 

t minCKs-l.is) 

^ j fi/,X I _ M _j + M + ,A' s -l,4^/,X,'_M_+f s ,i+M + -(A's-n+fs'l'' : s-f, 

,4=max(0,t s -l) 



P( (* - x n) 

(7b) 



because (x;_M_,/+M+\X;-M_+Ci+M + -(A: s -l)+4)u(x;_M_+4,(+M + -^ - (x;-M_,i+M + \ 

X/-m_+j(: s ,!+m + -a: s +*:,)■ Since (l6all6bl are valid for K s — 1 by Aitken's Lemma |Q1 pp. 204-205], (ITbl proves that they are 



also valid for K s = 2, and by induction VA^ e {l,--- ,M— 1 := M_+M + -l], completing the proof. Direct computation, 
using ( TTal i. proves that the consistency relation (l3bl holds for A" s = 1, and, by Lemma |2~T1 V/T s € {!,■•■ , M — 1}. 



Finally, strict positivity of B/, 



7.x,_ 



.1.1 



induction, using d6bb . V£ s e {0, ■ • • , K s ] and V^ s e {1, 



follows by the order relations assumed in (Eat , and then by 

,M~ 1 := M_ + M + - 1}. □ 

Because of the positivity of the numbers fi/.x,_„_ i+M ,A: s ,/t s £ R>o ( l6bl it is quite straightforward to study the sign 
of the weight-functions 07 Xj _ M _, +Jt ^^(jt) ( l6al >. which allows to determine the intervals on the real axis where the 
combination ((3) of the Lagrange interpolating polynomials on the substenicls is convex. 



Proposition 3.2 (Convexity in the neighbourhood of x,). Assume the conditions of Proposition l3.ll Furthermore as- 
sume that K s < |~y~|. Then the weight-functions of the combination (0) of the Lagrange interpolating polynomials on 
substencils (Proposition 13. It satisfy 



VM ± € Z : M := M_ + M+ > 2 
V* s e {!,•••, [f]} 



(8) 



V£ s e {0, 



Proof. Because of the consistency condition Obi , (non strict) positivity of the weight-functions o"/ iXi _ M j^jc^x) 
© suffices (proof by contradiction) to prove (|H). Rewrite (l6al l as 



A/. 



!-M_+/t s -l 



P[ (x- X„) 



n=i-M- 
i-M-+k s -l 



n 

i+M + 

n 



(x X w ) 



n=i+M+-K s +k s +\ 
i+M + 



\X Xfi) 



n=i+M+-K s +k s + l 



P| (x- X n ) 



k s =0 



0<k s <K s (9a) 



Obviously we have 

'i-M^+k s -\ 

sign Y\ (x - x„) 

i+M+ \ 

sign P[ (x - x„) 



= 1 Vk s 6 {1 



n=i-M- 



,K S ] Vx > max x,^ M + ^ s _i 

0<K<K„ 







V& s 6 {0, • ■ • ,K S - 1} Vx < min x i+M+ - Ki+ki+x = x, +A r + _s: s+ i 

0<L<K, 



(9b) 
(9c) 



Combining (|9al - l9cl l with the positivity of the numbers 6/, Xi 



! ->(l 



(16bt . and taking into account that 



(— l)*" s ^ s = 1, proves ((H). Notice that the condition for the interval [x,_m_+a;-i, ^i+M + -%ihi] in ((H) to contain at least 1 

cell (at least 2 grid-points) is -M_ + K s - 1 < M+ - ZT S + 1 <^^> 2K S < M + + M_ + 2 
which explains the additional constraint on K s included in the hypotheses of Proposition [3]2 



13 (A.2)] 



^s<|^| = |f | + 1, 
□ 



sJlSl). For cr ItX ._ ir _ n . +r ^ lA (x) the positivity interval is, by ©, [*;_(r-i)+(r-i)-i, x i+r _( r _i )+ i] 
[4, Tab. 2.1, 



Remark 3.3 (Typical stencils 

[Xj-uXi+i] in agreement with [4, Tab. 2.1, p. 507], while for cr !x . _ t ., r ,kS x ) the positivity interval is, by ((H), [x,_, +r _i, x, +r _ r+l i 
[x,_i, x, + i] in agreement with [4, Tab. 2.2, p. 511]. □ 



Lemma 12711 only requires the determination of weight-functions for the 1 -level subdivision. It is therefore not 
limited to a particular family of stencils and/or subdivisions, and can be used to determine weight-functions on biased 
stencils, eg near the boundaries of the computational domain. 
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Remark 3.4 (Relation to previous work). WENO interpolation applied to the development of central WENO schemes 
only requires the computation of the value of the weight-functions at specific points on the stencil, and these can be 
computed by solving a linear system [3]. Carlini et al. [2] pointed out that the weight-functions in representation (0) are 
of the form (l6al l. with unknown constants := (-l)*"' 4 '^^ ,+m + which can be determined, in the general case, 
by solution of a linear triangular system. These authors [2] studied in particular the stencil x,-_( r _iy + ,., for K s = r — 1, 
for which they obtained an analytical expression for the coefficients y^, and proved convexity Vjc e [xu x,+i]. Liu et 
al. J3l used the same form as Carlini et al. [2] for the weight-functions, computed the coefficients up to r — 6, and 
observed that the interval of convexity is actually [jc,-_i, (Remark |3. 3b . They also studied the stencil X,_ r i+r , for 
K s = r, gave an analytical expression for the coefficients y^, which were computed up to r — 5, and observed that 
the interval of convexity in this case is (Remark 13. 31 >. These results are in agreement with those proven 

in Proposition 13.21 (Remark l3~3t . Proposition 13.11 studies an arbitrary stencil X,_m_./+m + (Definition ll.il ). and level of 
subdivision K s e {1, • • ■ , M — 1} (M := M_ + M+), and obtains an analytical recursive expression for the coefficients 
fi/.x,_ M _,+M ,k s Jc m d6ab . In this way, we were able to give a formal proof for the interval of convexity, which was 
determined for general values of M+ and K s (Proposition l3.2l i. □ 

Remark 3.5 (Alternative formulation). The expression of the (K s = l)-level weight-functions (T7ab can also be written 
in an equivalent form, using ratios of fundamental functions of Lagrange interpolation on X,_m_,/+m + (Definition I l.U 
and on its {K s = l)-level Neville substencils (Definition ! 1.21 ). X,_M_,i+M + -i and X;_m_+i,;+m + - The Lagrange interpolat- 
ing polynomial of a function / : R — > M, on X,_ M can be expressed (ll (9.3,9.4), p. 184] as 

_ [1,(9.3)] \ "l , . „ N 

Pi,x,. M ^ i+u+ (x;f) = 2j ff 'A*.i*,i+fW/w ( 10a ) 



(=-M 
M + 

lil (9.4)1 l — r X — Xj+k 

:= I — V I e [-M-, ■ ■ ■ , M + \ (10b) 

J I Xi+i - x i+k 



k=-M. 



f M :=f{x i+e ) V t e [-M-, • • • , M + ) (10c) 



<Xl,X,-M 



where the M + 1 polynomials cti,x,_ M _ i+M ^,i+t 6 MmW dlObt are + 0r m [x](x), linearly independent, and form a basis of 
the space of all polynomials with real coefficients and degree < M, 1m[x] JSJ p. 2771]. It can be verified by direct 
computation, using definition ( II Obi ), that d7al ) is equivalent to 

, . UMMi «/,Xi_ M _, i+ M + ,i-M_(jc) 

,i,oW = rr (Ha) 

»/,X,-M_, + M + -l,!-M_W 



o"/,x,_„_, i+M+ ,i,iW = — (Hb) 



GilEob) Q '/,x,_M_,, + « + ,i+M + (-y) 

a /,x,_ M _ + , J+M+ ,i+M + W 

The interest of this alternative expression (Q~T), which is analogous to the expression of the {K s = l)-level weight- 
functions for the representation of the Lagrange reconstructing polynomial on a homogeneous grid [9, Lemma 4.2, p. 
2780], is that it can be generalized for the representation of the n-derivative of the Lagrange interpolating polynomial 
by combination of substencils, as will be shown in 2] □ 



4. Application to the n-derivative of the Lagrange intepolating polynomial 

One of the motivations that led to the formulation of Lemma l2.1l was the study of WENO reconstruction in view of 
the computation of f'(x), and this application is studied in [9] (results and relation to previous work are summarized 
in Remark |4~T1 >. The expression of the {K s = l)-level weight-functions for the Lagrange reconstructing polynomial 
is similar to (fTTT ). upon replacing the fundamental functions of Lagrange interpolation in dTTT > by the corresponding 
fundamental functions of Lagrange reconstruction ®, (32), Lemma 4.2, pp. 2780-2781]. It turns out that a similar 
relation is valid for the n-derivative of the Lagrange interpolating polynomial 

f=-AL 
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and these new results are formulated in Proposition 23] 



Remark 4.1 (Lagrange reconstructing polynomial). The case of reconstruction of a function h(x) from its cell-averages 
f(x), sampled on a given stencil (Definition ! 1 .It . is important for the construction of numerical schemes used in the 
solution of hyperbolic PDEs |2fl. In the particular case of homogeneous grids (x, + i -x, = Ax = const e R>o V/), recon- 
struction can be used for the computation of f'(x) fld] Q In the context of methods for the determination of numerical 
fluxes, what is needed are the values of the weight-functions at x i+ i (optimal weights [3]), which were usually com- 
puted from the solution of a linear system III (13), p. 8489]. Recently, Arandiga et al. lUOIl gave analytical expressions 
of the optimal weights for the (K s - r - l)-level subdivision of the upwind-biased stencil [i - (r — 1), ■ • • ,i + (r— 1)}. 
Liu et al. fl4J] showed that the weight-functions, in the reconstruction case, are rational functions, expressed the weight- 
functions for the (K s - \ 4H)-level subdivision of the stencil {i - |_4^J, ■ • • , i + M - |_-f J}, in the range M e {2, • • ■ ,11}, 
and computed the interval of convexity around x i+ i . In J9J], we use the recurrence of Lemma ETI to study in detail the 
weight-functions for the Lagrange reconstructing polynomial fl, obtain explicit recursive expressions for the weight- 
functions for an arbitrarily biased stencil on a homogeneous grid and an arbitrary subdivision level, and determine the 
interval of convexity in the neighbourhood of x i+ 1 . □ 

Liu et al. ]4J], also study the representation of the first two derivatives (n e {1,2}) of the Lagrange interpolating 
polynomial for the particular homogeneous stencils and subdivisions studied in the reconstruction case (Remark |4.U . 
In the present work, we show that the (K s = l)-level subdivision (Definition 11.2b weight-functions can be explic- 
itly determined for the n-derivative of the Lagrange interpolating polynomial (Q~2). Using Lemma |2~T1 we define the 
weight-functions for the representation of the n-derivative of the Lagrange interpolating polynomial (fT2t on X,_m_,/+m + 
(Definition ll.il ), by combination of the n-derivative of the Lagrange interpolating polynomials on the ^T s -level subs- 
tencils (Definition [Lj]), requiring that n < M - K s so that the n-derivative be + Or m k [ X ](x) on the substencils. The 
result is formulated in the following Proposition. 

Proposition 4.2 (Weight-functions for the n-derivative of the Lagrange interpolating polynomial). Assume the condi- 
tions of Proposition [XT] Then, VM ± e Z : M :- M_ + M + > 2, VK S < M - 1, Vn < M - K s , the n-derivative with 
respect to x of the Lagrange interpolating polynomial on X,_m_,i+m + (Definition ll.il ) can be represented, almost ev- 
erywhere, by combination of the n-derivative of the Lagrange interpolating polynomials on the K s -level substencils 
(Definition ! 1.21 ) o/x,_m_./+m+, as 

jfc s =0 

V/ : R — > K (13a) 
where the rational weight-functions 0"/ ()i)iX . _ M _ f+J/+ ,x s ,i s W are defined recursively by 



a (n) (x) 
a in) (x) 



K s = 1 



min(A" s -l,i s ) 

°" ;„ 11 ,x i _„_ iMt ,if,-u s (-t) o" ;„),x i _M_« s ,M + -K-i )t f„iA-4(- ( ) ^ s > 2 

( 4=max(0,ifc s -l) 

VJt s e {(),•■• ,K S ] \fK s 6{1,--- ,M-1} Vn e {0, ■ • ■ ,M-K S \ (13b) 
and satisfy the consistency condition 



I 

K=0 



(13c) 



on general inhomogeneous grids reconstruction does not provide f'(x), only numerical fluxes for the discretization of the PDE fjj 
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The set of poles of the rational weight-functions ,x,_ M _ /+M+ ,/r s ( 1 1 3at satisfies 

Swx,-„_„+„ t ,i < |vc R: afl M+u+M+ . +M+ (x) = o} 



K s -l L s 

L s =Of s =0 I L s =0 4=0 

WT S e {!,••■ ,M- 1} 



00=0 



(13d) 



(13e) 



Proof. Because of Lemma I2TT1 it suffices to prove that Proposition 14.21 is valid for K s = 1. Consider first the 
consistency relation ( 1 1 3cb for the (K s = l)-level weight-functions defined by d 1 3bb , 



a 



^/(„),Xj- M _, Wf+ ,l,oOO 



°~iM,Xi- M _ MM+ ,1,1W 



(n) 



(A') 



» 



g /,X,-M_,, +M+ ,'+M + ^ X ) 
a /,X;_M_ + u+ M + ,i+M + W 



Vn 6 {0, • • • , M - 1) 



Vn 6 {0, • • • , M - 1) 



(14a) 



(14b) 



By straightforward computation using the expression (1 10b) of the fundamental functions of Lagrange interpolation, 
we have 



M+-1 



Xi—M- ~ Xj+k 



( M+-1 ^ M+-1 



V<t=-M_ + 1 
f M+-1 



n 



X i+M+ ~ x i+k 



Lt±+ L 

n 



x Xj-t-k 

A:-l/ +1 X '' +M + ~ Xi + k 



a /,X i -M- + l,,'+M + ,'+'W + ( X ) 



(15a) 



U=-M_ + l 

and taking into account that the product in the last line of (II 5ab is independent of x, we have by differentiation 

( M+-1 



.(«) 

/,X,_m_j+M+ 



n 



« ( /"x V« e {0,--- ,M- 1} 



(15b) 



U=-M_ + l 

where « = implies no differentiation, ie (115al i. Since the product in (1 1 5bb is independent of «, we also have 



(n-l) / \ / \ 

g /.X,- M _., +M+ - 1 ,;-M- 1 -^ {Bb} ffi /,x,- M -+i, t M + ,i+M + W 

( n ) ( \ (A) ( \ 

Q; /,x,_„_.i + A f+ _i ,i-M- W a /,x,_„_ +1 ,, +M+ ,i+M + W 



Vn e {!,•■■ ,M- 1) 



(15c) 



We know by Aitken's Lemma [1, pp. 204-205] that the K s = 1 consistency relations hold for n — 0, and by direct 
computation from the expression (1 10b) of the fundamental functions of Lagrange interpolation, we have (ITal [TO 



a /,x i -„_,i +M+ ,i-*f-( J<: ) (MED *~*;+M + 



£*/,Xi-jf_,i + if + -i A-M- 00 



X i+M + Xi-M- 



«/,x,_j.,_,, +m+ ,!+M + (a) UMIB X-Xi-M_ 



°'/,X,-«_ +u+ „ + ,/+M + (a) 

which give by successive differentiation (proof by induction) 



(n) , ^ H5dl X - X,+m + j„) 



,_ m ,. v , J' - — ";" + _ M (x) - - 



n 



CP X ~ X i-M. („) . . 



■ffr„ ■ ... (^) 



(15d) 
(15e) 

(15f) 



..(n-l) 
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Using the above relations (fT4l[T5T l. we readily have 



0/<„>,x 



+ ,i,oW + °"/«,Xi_«_^ M) .,i,i(Jc) 



X i+M+ 



x X i-M_ 



X i+M + X i-M_ X i+M + X i-M_ 



= 1 



/ (k-1) 



X i+M + X i-M_ 



\n-i-) / \ (n— 1) / \ \ 

/,X,-„_, i+ „ + -, ,i-M_ W °7,X,-M- + l J+ M t W 



(«) 
\ a ix 



El, 



^0 



= 1 Vne {!,-•• ,M- 1) 



(16) 



By ( fTSI ). we have proven the consistency condition ( 1 1 3cb . for /f s = 1, Vn 6 {0, • ■ • , M — 1}, the case n — following 
from Proposition 13. II Obviously the weight-functions (TBI are defined almost everywhere, except at the roots of the 
denominator, which because of d!5bb . are defined by (II 3db . The c relation in (113dl > is used, because there may be 
cancelation of poles by polynomial division, as in the interpolation case (n = 0), where there are no singularities (|6}. 

To prove the representation (!13aK for K s = 1, we start from the remainder theorem of the Lagrange interpolating 
polynomial, which states [ 1, Theorem 9.2, p. 187] that for any real function / : K — > K of class C M+1 



V/eC M+1 (K) VxeK 3f(x;X,_ M _, /+M+ ;/)e [min(x, x,_ M J, max(x, x /+M+ )] : 

( M + A 



P/,x i _ M _, i+M+ fe/) =/(*) 



1 



(Af + 1)! 



j~[ (x-x /+i ) 



U=-M_ 



/( M + 1) W 



Since for any polynomial of degree < M we have 

Vq e EmW => | 9(4f+1)(jf) = Q 
we can write because of (I17ab 

M + -l 

P/,x i -M_ J+Jf+ _ 1 (*; g) =?W - ]~ [ (* - *>•+*) 



Vx 6 



k=-M_ 



Vtf e R M W : coeff[x M ,^(x)] = 1 



(17a) 

(17b) 

(17c) 
(17d) 
(17e) 



p/,x i _„_ +I , i+M+ (*; =?(*) - Pf (jc ~ Xi+k) 

i=-M_ + l 

Using the expression (llObb of the fundamental functions of Lagrange interpolation in (ll7clll7et PI and differentiating, 
we have, Vn e No, 

( M+-1 



(r, 9 )=« w (jc)- 



(«)/ 



=<? w (x) 



(x;q)=q w (x) 

r '.Xj_Jif_ + i,j + jf + v J/ J v ' 



(»)/ 



U=-M_ 



n * 

Vjfc=-M_ + 1 



a (n) fx) 



(17f) 



V^eRjvfW : coeff [x M , ? (x)] = 1 (17g) 



i-M- 



Xi+k) 



An) 



^,; + « + .i-M_( X ) 



(17h) 



M+-1 I M + -l 

k=-M- \k=-M- 



[ [ (.v - [ [ - x i+k ) I 07,x,-_ M _. /+M+ ,f-M_ (x) 



n 

U=-Jlf_+1 
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Combining ( [T4lll7fl[l7hb . and using ( I15bl ). we have 



^/<„>,x,.„_,, +M+ ,i,o(*) pfl_ M _ J+M ^ (x; q) + 0-/ w ,x I _ 4 ,_, i+M+ ,l,l(x) P?l_„_ tlM ^ (x; q) 



(o-/ ( „ ) ,x i _ M _,, +M+ ,i,o(^) + o- hnhXl _ M i+MtX i(x))q <n) (x) 



Bp, 

( M+-1 M+-1 

\\ (X i+ M + ~ X i+k ) P[ (x,-M_ - X i+ k) 



i=-M_ + l i=-M_ + l 



« r „ ■ i, (x) a,' ., (x) 







^ ^ +m (^ ;? ) v ? e R M [x] : coeff[x M , = 1 Vra e {1, • • • , M - 1} (17i) 
Using the representation (fTZt of the n-derivative of Langrange interpolating polynomials in dl7iK we readily have 

M + 

Z a /l-M_,, + « + ,^ (x) V< ? 6 Rm[x] : coeff[x M ,^(x)] = 1 Vn € {1, •• • ,M- 1} (17j) 

Applying ( |17j| i, successively, to the polynomials 

M + 

K«M 3 J~[ (x - x, +m ) = Vie {x,_ M _, • • • , x !+M+ ) \ {*,+*} Vfc e {-M_, ■ • • , +M+) (17k) 



m--M_ 
m±k 



yields 



a /?L-, i+M+ ,;-M>) = °W^) a ^x i -M-, i+ M + -i,i-M_W ^ <E3 ( 171 > 

^.^^W ^.X^^.l.oM ^L^^i+fW 

+^A JL «.uW <*£L_ +1 , i+4 , + ,^« W e !- M - + 1. • ' • . M + - 1} (17m) 



SB W».uW a £L_ +w , + W*) <=» G3B (17n) 

Replacing the expression ([T2l for the n-derivative of Langrange interpolating polynomials in (113 at . proves, by (11711 - 
H7nV (fl3at V/ : M — > M, for K s = 1. Lemma|XT]completes the proof. □ 



Remark 4.3 (Relation of Proposition l4.2l to previous work). To the author's knowledge, the case n > 3 has not been 
studied before. For n € {1,2}, the simple analytical recursive expression (113bt for the weight-functions agrees with 
the expressions obtained by Liu et al. [4], for the case of a homogeneous grid, using symbolic calculation for r e 
{2, ■ • ■ ,7}, for the (K s — r- l)-level subdivision of the stencils X,_, v+( ,_i) := {— r+l,--- , r— j} J3, Tab. 3.2, p. 516, for 
n = 1, and Tab. 3.8, p. 520, for n = 2] and for the (K s = r)-level subdivision of the stencils X,_ f , +r := {-r+i, ■ ■ • , r+ j) 
10, Tab. 3.5, p. 518, for n = 1, and Tab. 3.9, p. 521, for n = 2]. Proposition 14.21 formally proves the existence and 
analytical expression of the weight-functions, for the n-derivative (Vn e {0, ••■ ,M — K s }), for a general K^-level 
subdivision (K s e {1, ■ • ■ ,M — 1}), for an arbitrary stencil of M + 1 distinct ordered points (Definition II .1) . on an 
inhomogeneous grid. □ 
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5. Conclusions 



Every system of functions (l4at . depending on 2 integer parameters, which is equipped with an associated sys- 
tem of weight-functions satisfying a 1 -level subdivision property (I4bll4cl > also satisfies K s -level subdivision relations 
(Lemma [2.1l ). with weight-functions generated by the recurrence d4eb . which can be interpreted as an inverted gener- 
alized Neville algorithm pp. 207-208]. 

As an application of Lemma |2"T1 we developed simple explicit expressions for K s -level weight-functions of the 
Lagrange interpolating polynomial (Proposition 13. It on a general stencil in an inhomogeneous grid, which allow 
explicit determination of the interval of positivity of the weight-functions (Proposition 13. 2\ generalizing previous 
results Q[4]. By © the length of the positivity interval is M + - M_ - 2K S + 2 cells. 

We further investigated the existence of K s -le\el weight-functions for the n-derivative of the Lagrange interpolat- 
ing polynomial (n € {1, • • • , M - K^}). Having proved simple analytical expressions for the 1 -level weight-functions, 
Lemma 12.11 was applied to develop an analytical recursive expression for the K^-level weight-functions of the n- 
derivative (n e {1, • • • ,M - K s }). These results are valid for general inhomogeneous grids. 

Other potential applications of Lemma lXTI include WENO integration [4], the development of WENO schemes for 
biased near-boundary stencils, and other than Lagrange types of polynomial interpolation 
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